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Lattice results for the gluon propagator in SU (2) pure gauge theory obtained on large lattices are 
presented. Simulated annealing is used throughout to fix the Landau gauge. We concentrate on 
checks for Gribov copy effects and for scaling properties. Our findings are similar to the ones in 
the SU{3) case, supporting the decoupling-type infrared behaviour of the gluon propagator 
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1. Introduction 

Over the last years interesting results in Landau gauge gluodynamics have been found for 
gluon and ghost propagators (and, consequently, for the running coupling) both within the semi- 
analytical Dyson-Schwinger(DS) as well as Functional Renormalization Group (FRG) approaches 
[|l|, and with the help of lattice computations [|3[ ^ ||, ^, They have excited controversal 
debates about the behaviour in the deep infrared (IR) region. Within the DS and FRG approach it 
was demonstrated [§] that the behaviour strongly depends on the choice of the infrared limit of the 
ghost dressing function taken as a boundary condition for solving the (truncated) system of equa- 
tions. The so-called scaling solution exhibits an IR singular well-defined power-like behaviour of 
the ghost dressing function and correspondingly a vanishing gluon propagator in agreement with 
the quite attractive confinement scenarios invented some time ago by Gribov and Zwanziger on one 
hand and by Kugo and Ojima on the other. Moreover, it was in accordance with BRST invariance 
properties. 

Lattice results - as long as they are based on the assumption that one has to choose Gribov 
gauge copies as close as possible to the fundamental modular region - support with convincing 
numerical evidence the so-called decoupling solution with finite IR limits of both the gluon prop- 
agator and the ghost dressing function. Revisiting the case of SU (2) lattice gauge theory we give 
here further evidence for this observation by a consequent use of an efficient gauge fixing method, 
the simulated annealing algorithm. Since previous lattice investigations of the IR limit both in 
SU (2) as well as in SU (3) were carried out at quite strong bare coupling values in order to reach 
largest possible physical volumes the continuum limit was not really under control. Therefore, in 
the given contribution to LATTICE '09 we have a look into the scaling properties of the gluon 
propagator, which seem to be a bit more involved in the SU (2) than in the SU (3) case. We neglect 
so-called SU (2)-flips which enlarge the class of Landau gauge orbits and allow to extremize the 
gauge functional even further. 

2. Landau gauge fixing, Gribov ambiguity and gluon propagator 

In order to fix the gauge on the lattice we apply gauge transformations g{x) G G = SU (Nc), 
{Nc = 2,3) to sets {Ux^^} of link variables by mapping Ux.ji ^Ux,^. = g{x)Ux,^g\x + fi) . 
The set of all admissible {^Ux,^} for a given field {Ux,^} is called a gauge orbit. The Landau 
gauge d^A^ = for A^-+^/2,^ = {\/2iago){Ux.ji - t/vV)traceiess is fixed by searching for the local 
maxima gimx{x) of the gauge functional 

Fu[g] = ^Y,')^tTrmx,n. (2.1) 

In general for non-Abelian groups G more than one local maximum gimx{x) can be found, the 
so-called Gribov copies. Since the values of the gauge functional [g] and other gauge-variant 
quantities 0{^U) computed for various Gribov copies typically are correlated, further clarification 
of the gauge fixing condition is required. In case of Landau (or Coulomb) gauges for SU (Nc) gauge 
theories it was proposed ^ to choose the global maximum ggmxi^) among all local ones, thus 
introducing the fundamental modular region (FMR) inside the Gribov region, the latter defined 
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by the positivity region of the Faddeev-Popov operator [j8|]. In the present study we still keep the 
FMR condition. In practice it is hard - if not even impossible - to reach the FMR. Therefore, one 
is interested to improve the gauge fixing method and/or to apply the method of choice many times 
starting from random gauges in order to find the "best copies" closest to the FMR. 

In what follows we reconsider the question of how the gluon propagator depends on the Gri- 
bov ambiguity when the (Landau) gauge is fixed on a lattice. This question has been addressed 



for relatively small lattice volumes already in preceding publications (e.g. in |10|, |11|1). Note 
that in previous papers, instead of a comparison between two different gauge-fixing techniques, a 
comparison is made of the "best copy" (be) with respect to the maximal Fu [g] value achieved and 
"first copy" (fc), i.e. corresponding to a randomly chosen copy |T^, 12]. In [jlT]] "worst" copies 
were used for comparison, as well. There a visible Gribov copy effect was reported for the gluon 
propagator in the infrared. However, the effect appeared to be much more pronounced, when (i) 
the gauge orbit was extended by admitting nonperiodic (periodic up to Z{Nc)) gauge transforma- 
tions (flips) and when (ii) the standard overrelaxation method (OR) was replaced by a simulated 



annealing algorithm (SA) - always followed by finalizing OR steps [ ]13| , |14| , |15| ]. Moreover, the 
more efficient gauge fixing approach (SA -i- flips) led to a suppression of finite-volume effects, and 
indications were found that the influence of flip transformations on the gluon propagator gradually 
weakens with increasing linear lattice size L. But all these results were obtained at rather small 
lattice volumes. 

In the present study for SU (2) we neglect the flip gauge transformations and compare the SA 
method with the OR algorithm for strictly periodic gauge transformations. Our implementation of 
SA gauge fixing in the SU (2) case differs from that for SU (3) gluodynamics [^, ^ only in some 
technical details. We note that the "temperature" Tmax, from which SA cooling starts in the SU (2) 
case, is chosen to be the "critical" value Tcr of some phase transformation [[T^, ^, which takes 
place in the gauge-fixing field g{x) interacting with the Monte Carlo generated equilibrium field 



Ux,^ according to the gauge functional (2A). For the SU (2) case it looks like a higher order phase 
transition [ p^ ] or even like a "crossover". Anyhow, we have chosen T^ax = ^c/. = 1.1 in most cases. 
Simulated annealing, also known as "stochastic optimization method" ^VT\\, in principle allows 
getting arbitrarily close to the global maximum Fij[gg,„j(], if the number of SA "cooling steps" is 
large enough. This is the underlying idea of our "single copy" method successfully used in our 
SU (3) papers [||, In fact, in our SU (2) computations we have used long SA chains with O(IO^) 
steps from Tmax down to Tmin = 0.01 [^. 

The comparison for the unrenormalized gluon propagator obtained with SA versus OR tech- 
niques on 80"* lattices at jS = 2.30 is shown in Fig. [T], where we have plotted only results for mo- 
menta surviving the so-called cylinder cut [|l^, ^. One can clearly see a noticeable difference 
between SA and OR gluon propagator values in the deep infrared region <7^ < O.lGeV^, where the 
Gribov effect leads to a qualitative change of the behaviour of the gluon propagator. The question 
remains, whether the Gribov copy effect for the gluon propagator weakens with a further increase 
of the lattice volume. 

3. Check of scaling and multiplicative renormalizability 

The progress to reach the infrared regime had the price to consider lattice gauge fields on 
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Figure 1: Comparison of the unrenormalized gluon propagator obtained with SA and OR gauge fixing 
methods. The OR data are taken from Ref. [Btl. 
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Figure 2: The unrenormalized gluon propagator D{q^) for various (L, /3 ) pairs, i.e. at fixed physical volume. 



rather coarse lattices. To our knowledge the continuum limit expressed in a proper scaling and a 
multiplicative renormalization behaviour was not yet considered in detail on large volumes. We 
have made a step into this direction again neglecting the influence of Z(2) flip transformations. 

We have computed the gluon propagator on a sequence of lattices with increasing linear lattice 
size L and j8, choosing (L, j8) such that the physical volume was kept more or less constant. We 
have produced equilibrium ensembles of MC configurations and fixed the Landau gauge with the 
(single-copy) SA method for: (L,j8) = (40,2.2), (56,2.3), (80,2.4), (112,2.5), i.e. for a physical 
box size of approximately 10 fm. The bare gluon propagator D{q^) for these pairs of parameters is 
shown in Fig. ^ 

There are quite strong differences which have not been observed in the SU (3) case before. But 
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Figure 3: The renormalized gluon dressing function in the MOM scheme. 

this does not come unexpected. In the continuum or scaUng limit the bare propagator is expected 
to be multiplicative renormalizable as in perturbation theory. This means that between the results 
obtained at different lattice cutoffs a finite multiplicative renormalization up to lattice artifacts 
should be possible. In a finite volume - unavoidable for any lattice results - the multiplicative 
renormalization could be violated by finite-size effects. 

In accordance with the standard momentum-subtraction (MOM) renormalization scheme we 
have multiplicatively rescaled the bare gluon dressing functions Z{q^,p) = q^D{q^) (for all j8 val- 
ues considered) equating the renormalized values at some q = pLto the tree-level value ZrenilJ-^) = 1- 
The renormalization point was chosen at = 5.8 GeV^, sufficiently far away from the cutoff 
momentum = 79.4 GeV^ for j8 = 2.50. For illustration the finite renormalization factors 
Z(/i^,j8)/Z(/x^,j8 = 2.5) are shown in the Table |l[ They were obtained by interpolating between 
the 7 data points closest to the chosen scale /i for j8 = 2.5. 

Table 1: Finite renormalization factors for jx^ = 5.8GeV^. 





Z{^\I5)/Z{^\I5 = 2.5) 


2.2 


0.815 


2.3 


0.8925 


2.4 


0.9489 



The MOM-renormalized dressing function is plotted in Fig. ^. One clearly sees that the three 
curves for the renormalized dressing function Zrenil^) obtained at lattice sizes L = 56, 80, 1 12 lie 
nicely on top of each other, thus confirming the expected multiplicative renormalizability. For j8 = 
2.2 and L = 40 there are some scaling violations that can be understood as lattice artifacts. Some 
slight variations of the curve obtained with the largest lattice L = 112 probably can be attributed to 
problems with the still unsufficient Monte Carlo statistics and/or autocorrelations. 

4. Conclusions 

The comparison of OR- and SA-based results for the gluon propagator for L = 80 and jS = 2.3 
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clearly shows a noticeable Gribov copy effect in the range of momenta q'^ < O.lGeV^. At the 
moment the (dis)appearance of this effect for even larger volumes is an interesting open problem. 
An open question is also, whether our results obtained on large volumes will be modified, when 
Z(2) flips are taken into account. This is a matter of research in a forthcoming papers of one of the 



coauthors [19] 



Using SA-based Landau gauge fixing we have got numerical confirmation for a nonpertur- 
bative multiplicative renormalizability for the gluon dressing function. Finite-size effects already 
seem to be negligible in the given range of momenta. 
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